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CRITICAL  NUCLEATION  FIELD  AT  THE  STRUCTURED  SURFACE  OF  A  SUPERCONDUCTOR 


Z.  C.  Wu,  Daniel  A.  Jelski  and  Thomas  F.  George 
Departments  of  Physics  &  Astronomy  and  Chemistry 
239  Fronczak  Hall 

State  University  of  New  York  at  Buffalo 
Buffalo,  New  York  14260 


ABSTRACT 

The  effect  of  structured  surface  of  a  superconductor  on  the  critical 
nucleation  field  is  discussed  in  two  cases:  one  with  the  magnetic  field 
parallel  to  the  grating  wavenumber  and  the  other  parallel  to  the  ripples.  In 
the  first  case,  it  is  found  that  the  critical  field  is  reduced  as  a  function  of 
grating  height,  whereas  in  the  latter  case  it  is  increased. 
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1 .  Introduction 

The  purpose  of  this  paper  is  to  investigate  the  critical  nucleation  field 

(H  .,)  near  the  surface  of  a  Ginzburg-Landau  (GL)  superconductor.  We  assume  that 

c  j 

there  is  a  sinusoidal  grating  imposed  on  the  surface,  and  the  primary  purpose  of 
this  paper  is  to  investigate  the  effect  that  this  grating  has  on  the  critical 
field.  We  shall  consider  two  cases:  one  in  which  the  external  field  is  parallel 
to  the  grating,  and  the  second  where  the  field  is  perpendicular  to  the  grating. 

For  a  flat  surface,  this  is  an  old  problem,  having  been  previously  solved 
by  deGennes.^  The  existence  of  the  grating  at  the  surface  complicates  the  matter 
by  introducing  non-trivial  boundary  conditions,  and  we  shall  derive  a  solution  by 
solving  the  relevant  equation  (an  elliptic  equation)  numerically.  Thus  the 
primary  motivation  of  this  paper  is  to  extend  deGennes'  work  to  include  roughened 

surfaces.  A  secondary  motivation  is  to  discover  how  a  rough  surface  might  affect 

2 

the  properties  of  thin-film  superconductors,  especially  the  new  ceramic  variety. 
Unfortunately,  the  applicability  of  the  GL  theory  to  the  new  ceramics  is 
questionable.  The  coherence  length  appears  to  be  very  short,  and  since  the 
underlying  mechanism  is  not  known,  it  is  not  at  all  clear  whether  the  GL  theory 
is  the  proper  phenomenological  model.  Indeed,  recent  work  by  Chela-Flores  et  al 
constitutes  an  extension  of  the  GL  model  which  would  be  appropriate  if  the 
resonance  valence  bond  model  of  superconductivity  appears  to  be  accurate.  They 
claim  that  their  new  model  does  reproduce  the  observed  phenomena.  Another 
article  by  Lobb^  demonstrates  that  the  GL  theory  breaks  down  close  to  the 
critical  temperature  because  of  fluctuations. 

Therefore,  while  the  direct  application  of  the  current  work  to  the  new 
superconductors  is  problematical,  it  certainly  applies  to  traditional  materials. 
Further,  and  this  is  a  subject  of  ongoing  research  in  this  laboratory,  the 


methods  used  in  this  paper  are  probably  readily  extendable  to  the  work  of  Chela- 
Flores  et  al.  Finally,  given  that  this  is  an  obvious  extension  of  deGennes  work, 
we  consider  the  problem  to  be  of  intrinsic  interest. 

The  GL  theory  is  described  very  nicely  by  Tinkham,^  Schrieffer^  and  Parks, ^ 

S  9 

in  addition  to  the  original  paper.  Gor'kov  has  shown  that  the  GL  theory  can  be 

rigorously  derived  from  the  BCS  theory  under  the  condition  that  the  temperature 

is  close  to  T  ,  the  critical  temperature.  DeGennes  has  shown  that  under  the 
c 

appropriate  constraints,  to  be  described  later,  the  GL  equation  reduces  to  a 

Schrodinger  equation  for  a  harmonic  oscillator.  The  eigenvalues  of  the  equation 

can  be  solved  analytically,  and  hence  the  nucleation  field  in  the  bulk  (H^)  can 

be  determined.  In  the  reduced  units  that  we  are  using  (which  will  be  described 

in  due  course),  H  -  is  found  to  be  DeGennes^  has  shown  that,  for  a  flat 

c2  It r 

surface,  =  1/3.71.  We  shall  duplicate  that  result  in  what  follows,  and  also 
show  that  a  structured  surface  decreases  the  surface  nucleation  field  if  the 
magnetic  field  is  oriented  parallel  to  the  grating  wavenumber.  More 
interestingly,  we  will  also  demonstrate  that  the  critical  field  is  increased  for 
the  perpendicular  orientation. 

Section  2  is  a  description  of  the  GL  theory  theory  sufficient  to  define  the 
problem,  the  notation  and  the  units,  while  Section  3  contains  our  result  for  the 
case  where  the  external  field  is  parallel  to  the  grating  wavenumber.  Finally, 
Section  4  contains  the  most  interesting  result,  namely  the  case  where  the  field 
is  perpendicular  to  the  wavenumber. 


2.  Ginzburg-Landau  Theon 


In  the  Ginzburg-Landau  theory,  a  pseudo-wavefunction  <$(r)  is  introduced  as 

2 

a  complex-order  parameter.  Here  |$(r)|  represents  the  local  density  of 


superconducting  electrons,  n,,(r).  Near  the  critical  temperature,  the  free  energy 

2  2 

of  a  superconductor  can  be  expanded  in  powers  of  |$|  and  |V$|  .  The  condition 
that  the  free  energy  must  take  a  minimum  value  leads  to  the  Ginzburg -Landau 
equations. 


3|*|2*  +  2^(f  ?  -  “  l)1  *  -  0  , 


J  =  t-  curl  H  =  -  $7$*)  -  A  , 


(2.1) 


(2.2) 


where  m*  amd  e*  are  the  mass  and  charge  of  the  Cooper  electron  pair,  and  the 
parameters  (a  <  0)  and  0  are  temperature  dependent. 

A  derivation  of  the  Ginzburg-Landau  equation  by  the  variational  principle 
requires  certain  boundary  conditions  which  have  clear  physics  significance.  For 
an  insulating  surface,  the  boundary  condition  must  be1^ 


7  -  “  *)  *Lm  0  • 


(2.3) 


where  n  is  the  direction  normal  to  the  surface.  This  condition  assures  that 
there  is  no  supercurrent  in  the  normal  direction  of  the  boundary.  For  a  metal- 
superconductor  interface  with  no  current,  deGennes  has  given  the  condition^ 


<7  7  -  7"  ^  *L  =  7  *  • 


(2.4) 


where  b  is  a  constant. 


1/2  . 


If  <t  is  much  smaller  than  ,  where  4>  =  (-a/6)  is  the  value  of  the  wave 

OO  00 

function  in  the  deep  interior  of  the  bulk  material,  then  the  0-term  in  (2.1) 


5 


becomes  negligible  in  comparison  with  the  a-term.  For  this  case,  one  can  ignore 
the  6-term  and  the  Ginzburg-Landau  equation  (2.1)  is  of  the  form  of  a  linear 
elliptic  equation.  One  case  to  which  the  linearized  Ginzburg-Landau  equation  can 
be  applied  is  the  nucleation  of  superconductivity  for  bulk  material  with  or 
without  a  surface.  One  would  like  to  evaluate  the  highest  field  at  which 
superconductivity  can  nucleate  in  the  interior  of  a  large  sample,  or  near  the 
surface  of  a  bulk  sample,  in  a  decreasing  external  magnetic  field.  In  fact,  this 
is  the  best  case  for  linearization  to  be  valid.  Ideally,  at  the  critical  field 
the  superelectron  density  must  be  infinitesimal,  where  the  linearized  version  is 
accurate.  In  this  instance,  the  pseudo-wavef unction  must  decay  exponentially 
into  the  bulk.  The  nucleation  field  in  the  bulk  (H^)  and  near  a  flat  surface 
(H  _)  have  been  found  to  be^ 

CJ 


H 


c2 


tr)V 


(x|m  $ 

=  A-(  Ji') 

2/2' 


(2.5) 


and 


$ 

H  ,  =  1.695  H  „  =  1/1.371  (-4)  ,  (2.6) 

ci  cl  ^l 

Yic 

where  4>g  =  ^  is  the  quantum  flux  and  E,  is  the  Ginzburg-Landau  coherence  length. 

The  Ginzburg-Landau  theory  is  a  localized  theory  of  superconductivity, 
despite  the  fact  that  superconducting  phenomena  are  nonlocal.  In  the  Gor'kov 
derivation  of  the  GL  equation  from  BCS  theory,  the  nonlocality  is  represented  by 

the  symmetric  kernel  H  (r,r')»  connecting  two  vectors  r  in  coordinate  space, 

n 

where  the  frequency  =  (2n+l)T,  and  T  is  the  temperature  with  n  ranging  over 


m 


m 
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nonnegative  integers.  Then,  in  a  dirty  superconductor,  the  mean  free  path  of  the 
bulk  material  is  much  smaller  than  the  Pippard  coherence  length,  and  the 
coherence  length  for  the  symmetric  kernel  in  the  lowest- frequency  approximation 


_  _  f¥oa/2 

\xjQ  "  v6trT  ’ 


(2.7) 


where  V^,  is  the  Fermi  velocity.  On  the  other  hand,  the  Ginzburg -Landau  coherence 

1/2 

length  is  proportional  to  l/CT-T^)  ,  where  T^  is  the  critical  temperature. 

Therefore,  as  long  as  T  -*  T  ,  then  E,  will  be  much  greater  than  E,  and  we  can 

c  w0 

neglect  the  nonlocal  effect.  In  this  case,  the  deGennes  boundary  condition  (2.3) 
for  the  free  surface  would  be  valid  as  long  as  the  surface  structure  is  on  a 
scale  larger  than  E,  .  Since  the  Ginzburg-Landau  length  can  be  made  arbitrarily 

long  by  raising  the  temperature,  we  think  there  is  considerable  room  in  which  the 
results  derived  from  using  the  local  approximation  will  be  both  valid  and 
interesting,  and  it  is  this  case  that  we  consider  here.  In  general,  the  boundary 
condition  is  very  complicated  in  as  much  as  the  nonlocal  effects  have  to  be 


considered  explicitly. 


11,12 


Critical  Nucleation  Field  Parallel  to  Gratina  Wavenumber 


where  y  is  the  grating  amplitude.  For  the  first  case  the  external  field  H  is 
assumed  to  be  parallel  to  the  grating  wavenumber.  We  assume  z  to  be  the  distance 
from  the  structured  surface.  By  translation  symmetry  along  the  y-axis  the 
wavefunction  takes  the  form 

ik  v 

$  =  e  *  f(x,z)  .  (3.2) 


Then  the  linearized  Ginzburg-Landau  equation  becomes 


2itH  . 

(~7~ -)  (z-zn)  f(x,z) 
0 


f(x,z)  ,  (3.3) 

r 


where  we  use  the  gauge  A  =  (0,-H^z.O)  for  the  vector  potential  and  the  boundary 

condition  is 


9n‘ surface 


(3.4) 


In  order  to  obtain  the  maximum  value  of  H  ,  for  which  f(x,z)  has  a  bound 

x 

solution,  one  has  to  solve  the  elliptic  equation  numerically.  Then  zq  is  an 
adjustable  parameter,  and  it  related  to  the  y-component  of  the  supercurrent 
expressed  a S 


=  (tS  k  +  zJt  H  x) 


(3.5) 


we  solve  the  elliptic  equation  in  the  over  one  period  of  the  grating  as  shown  in 

13  14 

Fig.  1.  The  differential  equation  is  replaced  by  a  difference  equation.  ’ 


Assuming  periodic  boundary  conditions,  the  value  of  the  wavefunction  at  boundary 
AO  is  identical  to  that  of  boundary  CD.  At  the  surface  (boundary  AC)  we  use  the 
von  Neumann  conditions,  i.e.,  the  normal  derivative  of  the  wavefunction  must  be 
zero.  We  assume  that  boundary  OD  is  sufficiently  far  enough  away  from  the 
surface  so  that  surface  effects  can  be  ignored.  We  can  thus  use  the  analytic 
solution  for  the  flat  surface  as  the  boundary  condition  here;  the  value  of  the 
wavefunction  at  this  boundary  can  be  an  arbitrary  constant  (normalization 
condition) . 

We  arrive  at  the  maximum  value  of  H  for  the  linearized  GL  equation  by 

x 

testing  the  convergence  at  the  bottom  boundary  (deep  in  the  interior),  and  also 

the  behavior  near  the  surface.  This  result  corresponds  to  the  bound  state  of  the 

minimum  excitation  for  a  given  z^.  However,  z^  is  an  adjustable  parameter,  and 

one  has  to  vary  zn  in  order  to  maximize  the  resulting  value  of  H  .  In  practice, 
u  x 

this  results  in  a  two-dimensional  variational  problem,  with  the  adjustable 
parameters  being  Hx  and  Zq.  We  find  that,  especially  for  large  gratings,  the 
solution  can  be  very  unstable  to  small  variations  in  these  parameters,  and  hence 
the  calculations  become  quite  time  consuming.  Fig.  2  illustrates  a  typical 
result. 

Figure  3  shows  the  numerical  result.  We  have  obtained  the  nucleation 

critical  field  for  the  grating  wavenumber  about  10?  ^  and  a  grating  amplitude 

ranging  from  0  to  1.2  GL  coherence  lengths.  When  the  grating  amplitude  is  zero, 

1.70  $ 


we  find  =  l-^OH^ 


-,  which  is  consistent  with  the  result  obtained  by 


2tt5 


de  Gennes  for  a  flat  surface.  As  the  grating  height  increases,  decreases  and 


asymptotically  tends  toward  the  value  for  the  bulk  material,  H  7  =  - x-  This  is 

cz  2ir  £ 

expected  since  as  the  grating  height  tends  to  infinity,  the  structure  of  the 
surface  becomes  a  series  of  superconducting  thin  films,  oriented  now  parallel  to 


r  K."  X.-  v~  ht.  )i.t  nr.  )C1 


the  z-axis.  The  direction  of  the  magnetic  field  is  in  the  x-direction,  and  hence 
is  oriented  vertically  with  respect  to  the  thin-film  layers.  The  wavefunction 
must  be  the  bound  state  with  minimum  excitation  which  satisfies  (3.5).  If  k  =  0, 
then  the  solution  satisfies  the  von  Neumann  conditions  for  the  thin  film.  Thus 
the  critical  field  must  be  equated  to  that  for  the  bulk  material,  •  This 
ignores  the  effect  in  the  z-direction,  which,  however,  only  affects  the  region 
very  close  to  the  "top"  of  the  grating.  Since  the  grating  amplitude  is  infinite, 
this  effect  can  be  ignored. 

It  is  also  apparent  that  the  critical  field  does  not  change  significantly 
for  A  <  0.3£.  This  is  because  the  von  Neumann  condition  forces  the  wavefunction 
to  be  constant  near  the  surface.  If  the  amplitude  is  much  smaller  than  C>  then 
this  condition  implies  that  the  solution  to  the  GL  equation  can  not  vary 
significantly.  As  the  grating  becomes  deeper,  however,  on  the  order  of  £,  then 
the  ripple  forces  a  dramatic  change  in  the  value  of  the  wavefunction. 

4,  Critical  Nucleation  Field  Perpendicular  to  the  Grating  Wavenumber 

Another  interesting  case  we  now  consider  is  the  nucleation  field  in  the 
ripple  direction  ( y-direction) .  The  gauge  for  the  vector  potential  is 


1  =  (Hyz,  0,  0)  . 


(4.1) 


Then  the  Ginzburg-Landau  equation  can  be  written  as 


82f(x,z)  _  8^f (x,z)  3f(x,z) 

8x2  8z2  *0  8x 


2irH  z  _  Tin 

+  ( — T*—)  f(x,z)  =  ( - -k2)  f(x,z)  , 

*0  ir  y 


(4.2) 


TO! 


where  we  have  used  the  ansatz  (3.2)  due  to  translational  symmetry  along  the  y- 
axis.  One  can  rescale  the  factor  (-  -  k“)  of  the  right-hand  side  of  (4.2) 

r  y 

2itHv 

to  1 ,  in  which  case  the  quantity  -  —  -  on  the  left-hand  side  becomes  = 

o 

2itH  / <t>  (IsJbd.  _  k2).  if  one  finds  the  maximum  value  H  at  which  a  bound 
y  o  ^2  y  yc 

solution  of  (4.2)  exists,  then  the  corresponding  critical  nucleation  field  is 


(4.3) 


with 


k  =  0,  (4*4) 

y 

which  is  similar  to  the  previous  cases. 

We  numerically  solve  the  elliptic  equation  (4.2)  by  the  same  technique 
described  in  Sec.  3.  The  only  difference  is  the  boundary  condition  at  the 
structured  surface,  written  as 


3_f 

3n 


ie 


w  A  f 
nc  n 


(4.5) 


Figure  4  shows  the  numerical  result.  We  have  obtained  the  nucleation  crtical 
field  H  ^  for  the  grating  wavenumber  about  10£  *  and  a  grating  amplitude  ranging 
from  zero  to  0.6  GL  coherence  lengths.  Of  course,  when  the  grating  amplitude  is 
zero,  we  recover  the  deGennes  result  for  a  flat  surface.  Unlike  the  case  for 


which  the  external  field  is  parallel  to  the  grating  wavenumber,  here  the 
calculation  shows  that  the  critical  field  increases  as  a  function  of  grating 


anqnaronowgKMvi.>  ■.v.vva’v.v  s.’  -.  ■  =a  yr  >m«  v  k.«  e»  ^■owr;r^.nw,v.xv.?^^.v  - 
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amplitude.  It  is  expected  that  the  value  of  the  critical  field  will  reach  the 
value  for  superconducting  films,  with  a  thickness  of  the  same  order  as  the  period 
of  the  structure.  It  is  noted  that  if  the  thickness  of  a  film  is  much  smaller 
than  the  London  coherence  length,  the  critical  field  must  be  much  greater  than 
the  bulk  value. 
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Figure  Captions 

Fig.  1  A  schematic  drawing  of  the  system  under  discussion.  Letter  0  refers  to 
the  origin  of  our  coordinate  system.  The  nature  of  the  boundaries  is 
described  in  the  text.  The  amplitude  of  the  grating  is  y,  whereas  the 
wavelength  is  0.65. 

Fig.  2  This  shows  the  wavefunction  at  the  nucleation  critical  field  for  a 

2it  0.26*q 

structured  surface.  Here  y  =  0.225,  v  =  o~65  anc*  ^c3  =  - 2 —  an<*  zq  = 

0.645.  While  the  amplitude  of  the  wavefunction  is  determined  by  an 
arbitrary  normalization  constant,  close  inspection  shows  that  the 
magnitude  of  the  electron  density  decreases  in  the  ripples. 

Fig.  3  Nucleation  critical  field  vs.  grating  amplitude  for  a  field  oriented 

parallel  to  the  grating  wavenumer.  The  amplitude  is  measured  in  units  uf 

<t> 

0 

the  GL  coherence  length  whereas  the  critical  field  is  in  units  — .  The 

wavenumber  of  the  grating  is  about  10/5  (see  text). 

Fig.  4  Nucleation  critical  field  vs.  grating  amplitude  for  a  field  oriented 
parallel  to  the  ripples.  Units  are  the  same  as  in  Fig.  3. 
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